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A survey is presented of recent measurements in jet physics, and improved determinations of the QCD coupling 
constant Qs that these have made possible. 



1. INTRODUCTION 

Following the first realisation that hadronic 
particles consist of charged fermions known as 
quarks, held together by bosonic particles known 
as gluons, it was learnt that the quarks and glu- 
ons ("partons") can never be observed on their 
own. Any attempt to isolate a quark, for exam- 
ple by knocking it violently out of a hadron in an 
energetic collision process, does not result in an 
observable quark (or gluon). Instead, a shower 
of elementary particles is seen, referred to as a 
jet, containing the energy of the original quark 
or gluon. The gluons couple to the quarks by a 
running coupling constant as- 

In this account, we first make some brief notes 
regarding QCD processes. There follows a sur- 
vey of recent experimental results regarding jet 
measurements and what has been learned from 
them. In particular, we then examine some re- 
cent measurements of as and attempt to make 
some comparisons. 

2. STATUS OF CALCULATIONS 

The most basic diagram in deep inelastic scat- 
tering (DIS) consists of an incoming lepton ex- 
changing a virtual boson (photon or electroweak) 
with a quark in an incoming hadron (Fig. [l](a)). 
The quark is ejected to give rise to an observ- 
able jet. The coupling of the virtual boson to 
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Figure 1. Examples of hard QCD processes at 
order zero and unity in as- 



the lepton and the parton is electromagnetic or 
electroweak, and as is not directly involved. 

Higher order processes involve the radiation of 
gluons at various points. The outgoing quark may 
radiate a gluon, or an incoming quark may do so, 
in which case the exchanged boson may interact 
with it by means of a quark pair, a process known 
as boson gluon fusion. These processes are illus- 
trated in Figure [Ijb), (c) and are of order unity 
in as- Processes of higher order in as have also 
been studied and will be discussed below. 

At present fT\, all the relevant QCD processes 
at lowest order (LO) in as have been calculated, 
and also all the next-to-leading order (NLO) pro- 
cesses involving two partons entering and two 
leaving. Most of the NLO processes involving 
the radiation of a third particle have also been 
calculated, but little has been done at higher or- 
der. The NLO calculations, as we shall see, are in 
good agreement with the relevant measurements 
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at HERA. However theory errors are beginning to 
become dominant in the ag determinations, rais- 
ing the issue of how to obtain increased theoreti- 
cal accuracy. 

3. EXPERIMENTAL METHODS 

The experimental methodology for measuring 
jet cross sections is by now, in principle, well es- 
tablished. It is taken that the jet corresponding 
to a quark or gluon radiated in a hard process, 
namely at high transverse momentum px, will 
comprise particles that are in a definable sense 
close together. Two main approaches are now 
adopted: to define the closeness in terms of a 
cone angle in 77, space, where 77 is pseudorapid- 
ity and 4> is angle, or to use a clustering method. 
With cones, a cone radius is specified and par- 
ticles are accepted if they are all found in one 
cone. Subtleties arise in connection with the ba- 
sic cone search, and in cases where jets overlap. 
With clustering techniques a parameter is defined 
involving the transverse energies of two objects 
and their angular separation in r/, 0, and cluster- 
ing is performed iteratively between objects (in- 
cluding already clustered objects) whose cluster- 
ing parameter is less than a given value. 

The complicating issues here arise in connec- 
tion with the comparison with theory. In theo- 
retical calculations at a given order we will just 
have distributions of partons, and it will be de- 
sired to apply the same jet algorithm to these to 
determine the jet properties at the parton level. 
If only one parton were to be involved there is no 
problem. However this will not be the case owing 
to the radiation of extra partons, and these pro- 
cesses tend to have infra-red and coUinear diver- 
gences; the jet algorithm must not be disturbed 
by these if we are to have a definable correspon- 
dence between the parton jets and the hadronic 
jets, with just a hadronisation correction to be 
evaluated and applied. 

The effects in question tend to be small in mag- 
nitude and were ignored for many years in the 
context of cone jets, but this is no longer ac- 
ceptable if we are dealing with high order pro- 
cesses or precision measurements. Clustering ap- 
proaches are much less problematic, and have be- 



come standard practice at HERA; however with 
many-jet events the well-defined radius of a cone 
can give advantages. The mid-point algorithm, 
which performs a second iteration in its cone 
search, using seeds midway between the jets al- 
ready found, is helpful but not perfect. Recent 
developments here include SISCone (Seedless In- 
frared Safe cone) and the so-called "anti-kT" al- 
gorithm. These approaches should serve well at 
the LHC. For a more detailed account of these is- 
sues, the reader should consult the account here 
by C. Soyez [2]. 

The experimental approach, then, is to use a 
jet finder that will give an adequate relationship 
between experiment and theory, and check that 
the basic kinematic distributions of jet quanti- 
ties agree with the QCD calculation of choice. It 
will also be necessary to use a trusted set of par- 
ton densities (PDFs) in the proton. Discrepancies 
may indicate the need for more accurate PDFs. If 
all is well, the data may be used to extract a value 
of as, provided that the process studied involves 
diagrams that are of at least order unity in a^. 
A high-order process would give greater sensitiv- 
ity but may well not have been calculated. These 
tasks, of course, may not be completely separable. 

4. TEVATRON RESULTS 

The highest jet energies at present are achieved 
at the Tevatron operating pp collisions at a cen- 
tre of mass energy of 1.96 TeV. Both CDF and 
DO have made extensive studies in the field of 
jets and their properties, and we show here some 
of their recent results [3|4j . Figure [2] shows in- 
clusive jet distributions as function of transverse 
energy for different rapidity ranges. Comparison 
is made with an NLO QCD calculation (NLO- 
jet) and excellent agreement is observed, both in 
the shape of the distributions and in the absolute 
magnitude, in all ranges. The agreement in mag- 
nitude amounts largely to a confirmation of the 
CTEQ PDFs used, while the agreement in Et 
shape points to the accuracy of the QCD calcu- 
lations. 

Both collaborations have also presented cross 
sections for prompt photon production, which 
tests different perspectives on QCD. CDF find 
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Figure 2. Jet distributions as a function of trans- 
verse energy from CDF and DO, using the mid- 
point algorithm with cone radius 0.7. 



good agreement with theory ( JETPHOX) in their 
results, while DO find some discrepancies with 
theory in their photon plus jet cross sections. The 
cause for these is not at present clear. CDF's di- 
jet mass spectrum, measured up to 1400 GeV/c^, 
again shows excellent agreement with NLO the- 
ory (NLOjet + CTEQ6.1M) with no signs of any 
additional structural features. 




e / Y 




Figure 3. Direct and resolved processes at LO in 
photoproduction. 



5. JETS IN PHOTOPRODUCTION 

Figure [3] indicates the two main categories of 
process in hard photoproduction, at lowest or- 
der in as- In direct processes the photon cou- 
ples directly to a high-p^ quark line, while in 
resolved processes, the photon interacts through 
a hadronic intermediate state to which it cou- 
ples non-perturbatively; this state represents the 
"hadronic structure of the photon" and the hard 
scatter, describable perturbatively in QCD, takes 
place with one of the partons in this structure. 
The distinction between direct and resolved pro- 
cesses can be made also at higher order, but an 
arbitrary division must be made between them; 
only quark pairs with px above the so-called fac- 
torisation scale will be considered hard, while 
quarks at lower px values will be considered as 
constituting part of the hadronic structure of the 
photon. 

The fraction x-y of the photon energy that is 
taken by the jet system, for two or more jets, can 
be calculated from the kinematics of the jets. For 
a direct process it should ideally take the value of 
unity, but owing to hadronic effects and soft radi- 
ation its observed value is less than this. Direct- 
dominated and resolved-dominated processes can 
largely be separated experimentally by taking Xj 
above and below a value of approximately 0.75. 

ZEUS have published improved measurements 
of dijet production at HERA [5]. The shape of 
the Et distribution is well described over more 
than three orders of magnitude (Fig. 2]), while 



4 



P. J. Bussey 




Figure 4. Distribution of the mean Et of 
photoproduced dijets for direct-dominated and 
resolved-dominated event samples. 



Figure 5. Azimuthal angle between the jets in 
the analysis of Figure 21 compared with different 
theoretical approaches, for direct-dominated and 
resolved-dominated event samples. 



the predictions using two different hadronic mod- 
els of the photon |6|7j differ slightly for the re- 
solved events; of course there is little difference 
between the photon models for the direct events, 
where the excellent agreement between theory 
and experiment confirms both the general the- 
ory and the choice of proton PDFs (CTEQ). The 
azimuthal angle between the jets, which at LO 
should be 180°, is sensitive to summed higher or- 
der effects including soft radiation. The hadro- 
nising Monte-Carlo HERWIG is found to model 
this distribution better than the standard param- 
eters of PYTHIA, although the latter can in fact 
be tuned. (Fig. ^ 

ZEUS have also studied multi-jet distributions 
with a view to checking the accuracy of higher or- 
der calculations, and also to ascertain whether a 
multi-parton description of the resolved photon- 
proton scatter should be employed. The distri- 
butions in Et of the jets in three-jet and four-jet 
events, in various rj intervals, are found to be sat- 



isfactory. There are some indications that mul- 
tiple parton interactions (MPl) improve the de- 
scription. However clearer indications come from 
the x-y distributions (Fig. For high masses of 
the jet system, direct processes dominate, but for 
low masses there is a substantial resolved compo- 
nent which is not well described without the MPI 
option to PYTHIA or HERWIG. Even here the 
agreement is not perfect, but at least the magni- 
tude is well given. The effect is more pronounced 
in the four jet events, but is evident in the three- 
jet sample. 

As a further study in photoproduction, ZEUS 
have used the angular correlations in three-jet 
events to compare the predictions of standard 
QCD with those from various alternative models 
[8]. The following parameters can be defined: 

• 6h — angle between plane of highest energy 
jet and beam, and plane of two lowest en- 
ergy jets. 

• a23 — angle between two lowest energy jets. 
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Figure 6. Distributions in the observed value of 
x-y for three- and four-jet events for jet masses 
below and above 50 GeV. 



• Pksw = Korner-Schierholz-Willrodt angle 
(based on cross products of jet vectors). 

• '^max — maximum pseudorapidity of the 
three jets. 

Distributions of these variables can then be com- 
pared with the predictions from different variant 
models. QCD makes use of the symmetry group 
SU(3); we can compare this with SU(N) with 
large N, U(l)'^ (no triple gluon coupling), and 
SO (3). The results, unsurprisingly but gratify- 
ingly, favour standard QCD. Figure [7] illustrates 
some of the distributions. In each case, standard 
QCD fits the data well, while the other possibili- 
ties all fail in some regions of some of the plots. 

6. EFFECT OF POLARISED BEAMS 

Turning to deep inelastic scattering, we note 
that HERA II was able to deliver longitudinally 
polarised electron and positron beams. In the 



Figure 7. Distribution of angular parameters of 
the three-jet system with predictions from stan- 
dard QCD (solid line) and several variant theo- 
ries. 



electroweak sector of the standard model, the 
charged-current (CC) jet cross sections are sen- 
sitive both to the polarisation and to the nature 
of the lepton beam. Since CC DIS converts the 
incoming lepton into a neutrino, the signature for 
this process is at least one pT-imbalanced hard jet 
in the detector. The cross sections become mea- 
surable only at high values of the virtual boson 
mass squared Q^. As in all the HERA measure- 
ments, a kx cluster algorithm was used for jet 
identification. 

Figure [5] illustrates the results obtained [S] • It 
is clear that the cross sections are much higher 
for positive polarisations than negative for a 
positron beam, and vice versa for an electron 
beam. Agreement with theory, using a ZEUS 
parameterisation of the proton PDFs, is very 
good, and was confirmed using the CTEQ6 and 
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ZEUS CC e*p 
• P=-U.37 176 pb 'l 
o P=+0.32(lll2pb" 
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Normalised Inclusive Jet Cross Section 



Figure 8. Cross section measured by ZEUS for 
inclusive jets in charged-current DIS for positron 
and electron beams, and opposite longitudinal 
electron/positron polarisations, compared to SM 
theory at NLO (MEPJET) as a function of jet 
pseudorapidity. 
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Figure 9. Cross section for inclusive jets in the 
Breit frame, with Et > 7 GeV. The values are 
normalised to the total neutral current cross sec- 
tion at a given value of Q^. 



MRST2001 PDF sets. Comparable results are ob- 
tained for the cross sections as a function of Et of 
the jets, and for the cross sections evaluated for 
zero beam polarisation. Here we are specifically 
testing electroweak theory in a new regime from 
the measurements at LEP. 

7. BREIT FRAME MEASUREMENTS. 

A good method for measuring as in DIS is to 
determine the cross sections for different kine- 
matic ranges as a function of and make use 
of the DGLAP equations, which involve as, to 
describe the variation of the cross sections with 
and the kinematic variables. A fit to the entire 
data set enables as to be extracted. 

Here, however, we concentrate on ag determi- 
nations which depend more explicitly on measur- 
ing jet cross sections. HI have performed this 
task as follows. First, they take each event as 
viewed in the Breit frame and find jets using the 
kx cluster method in this frame. In the Breit 
frame, the event axis is taken along the direc- 
tion of the virtual photon, such that in a simple 
deep inelastic scatter (Fig. 1(a)) the quark is hit 
head-on and merely reverses its direction. Such 
scatters, although they are hard, have no trans- 



verse energy in the Breit frame. It follows that 
high-i?T jets in this frame result exclusively from 
higher order processes, namely those that involve 
non-zero powers of as. The sensitivity to as is 
therefore enhanced by measuring jets in this way. 

A laboratory rapidity cut is imposed in order 
to have well-defined acceptance for the jets, and 
a number of different differential cross sections 
are evaluated. Figure [H] shows normalised inclu- 
sive jet cross sections measured under these con- 
ditions, compared to a NLO calculation; similar 
distributions were obtained for two- and three-jet 
final states. The normalisation approach enables 
some of the experimental and theoretical uncer- 
tainties in as to be reduced. This analysis was 
performed at high values of Q^, which reduces 
some of the theoretical errors at the cost of poorer 
statistics. All the control plots gave good fits us- 
ing the NLO theory. From a general fit to the 
data, as was evaluated as a function of Q^. 

In a second analysis, HI used a lower range of 
Q^, giving better statistics but, as we shall see, 
notably worse theoretical uncertainties. Events 
are selected in the range 5 < < 100 GeV^ with 
the fractional virtual photon energy in the range 
0.2 < y < 0.7. Jets are selected, as before, in the 
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Breit frame and must have Ex > 5 GeV. Detector 
effects are simulated and acceptances evaluated, 
as before, using standard Monte Carlos. General 
comparison is made with the NLOJET++ Monte 
Carlo predictions. 

As is evident in Figure llOi there is a notice- 
able difference between the predictions for dif- 
ferent choices of renormalisation scale. Both 
and {E^ + (5^)/4 are plausible, but the former 
value provides a good fit to the data while the 
latter is barely acceptable. The difference be- 
tween the final results obtained with these two 
choices must therefore be accounted as a theoret- 
ical uncertainty. Varying the factorisation scale 
likewise gives rise to a theoretical uncertainty on 
the results. The experimental systematics must 
likewise be included in the normal way, and vari- 
ations in the chosen PDFs must also be taken. 

For a series of bands, double differential 
cross sections (Pa/dQ^dEr were evaluated for a 



set of jet Ex values. Fits varying Ofg were carried 
out on these. The results are given in the next 
section. 

8. RESULTS FOR 

Experimental values for as can be determined 
in various ways, and they are seen to fall with 
increasing values of the typical momentum scale, 
such as Q or an equivalent variable, demonstrat- 
ing the running nature of this coupling constant. 
It is conventionally quoted at = M|. At LEP, 
some accurate determinations have been made us- 
ing the properties of event shapes. Dissertori et 
al TI)] have used an NNLO QCD calcination to 
fit ALEPH data and obtain 



a, 



(Mz) 



0.I240±0.0008±0.00I0±0.00II±0.0029 



where the errors are, in order, statistical, sys- 
tematic, hadronisation and theoretical. From an 
analysis of ALEPH and OPAL thrust data, and 
incorporating their own theoretical model, Becher 
and Schwartz [llj have obtained 
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Figure 10. HI differential cross section values for 
inclusive jets in the Breit frame in DIS for a lower 
range of values. Left: standard renormalisa- 
tion scale; right: variant renormalisation scale. 



a, 



s(Mz) 



0.1172±0.0010±0.0008±0.0012±0.0012 



Using an NNLO -|- NLLA ansatz, Bethke et al. 
[n\ obtain from the JADE data 



= 0.1172±0.0006±0.0020±0.0035±0.0030. 



These values have high numerical precision but 
have potential dependencies on the theoretical 
approach that is adopted. 

From HERA we first examine the recent HI 
high-Q^ results, obtained from the data described 
above (Fig. Ilip , and representing a clear improve- 
ment on their publication of 2007 [13] . Sensitivity 
to Us is increased by measuring events with dif- 
fering numbers of jets. Having demonstrated that 
Us runs as expected with Q^, the group perform 
a combined fit and obtain the value 



^s(Mz 



0.1182 ± 0.0008 tomti ± 0.0018. 



The low-Q^ results give a fitted value of 
= 0.1186 ±0.0014 



^s(Mz 



0.0132 I n r,Q21 
0.0101 ^ u-uuzi 



where the three uncertainties correspond to ex- 
periment (total), theory, and PDFs. These results 
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Figure 11. Fitted values of from HI data 
at high Q^, as a function of Q. The narrow 
central band indicates experimental uncertainties 
and the broader outer band indicates theoretical 
uncertainties. 



(Fig. [T2I top) have a very much larger theoretical 
uncertainty than the high-Q^ set. However when 
the values obtained are plotted with the fitted 
values and certainties obtained from the high-Q^ 
set extended down to lower Q^, the agreement 
is perfect, indicating the consistency of the cen- 
tral values of the theoretical parameters applied 
over the whole range of Q^. Given the smallness 
of the experimental compared to the theoretical 
uncertainties, the need for NNLO calculations is 
becoming increasingly evident. 

ZEUS have evaluated as from some recent 
measurement of inclusive jets in photoproduction 
(Fig. [HI), taking the MRST2001 and GRV-HO 
proton and photon structures as central values. 
They obtain 



^s(Mz) 



0.1223 ± 0.0001 ± 0.0022 ± 0.0030, 



where the uncertainties are respectively statisti- 
cal, systematic and theoretical. This is compet- 
itive with other values and in good agreement 
with them. In DIS, HI and ZEUS have put to- 
gether their results for HERA I inclusive jet data 
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Figure 12. Fitted values of as from HI data 
at high Q^, as a function of Q. The narrow 
central band indicates experimental uncertainties 
and the broader outer band indicates theoretical 
uncertainties. In the lower plot the high-(5 band 
is extrapolated into the low-Q region. 



to produce a combined fit for as- Apart from 
the improved statistical accuracy, this approach 
allows for a correct account to be taken of sys- 
tematic and theoretical effects which are com- 
mon or correlated between the two experiments. 
From HI, 24 data points at different Et are taken 
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Figure 13. Inclusive jet distribution measured by 
ZEUS in photoproduction with Et > 17 GeV, 
compared to NLO theory. 



with 150 < Q2 < 15000 GeV2, while ZEUS con- 
tributes six inclusive cross sections in the range 
125 < Q2 < 100000 GeV2. The QCD cross 
sections are calculated with NLOJET++. The 
factorisation scale is taken as Q, the renormali- 
sation scale is taken as Et of the jets, and the 
MRST2001 PDF sets are employed. A matrix 
is defined and minimised using a Hessian method. 
The result is: 



s(Mz) 



0.1198 ± 0.0019(ea;p) ± 0.0026{theory), 



competitive with the LEP results. 

Figure fT4l illustrates results from the combined 
fit, showing the running of with jet Et, and 
also a collected set of ctg values from a number of 
HERA analyses, again showing well the running 
of as with the relevant momentum scale. A collec- 
tion of determinations of as is shown in Figure [TSl 
in comparison with a recent world average value. 
It is clear that the different approaches are in 
general in good agreement. The tendency is also 
evident that in order to obtain small theoretical 
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Figure 14. Fitted values of ag from HI and ZEUS 
data at high Q^, as a function of Q. The top plot 
illustrates results from the combined fit; the lower 
plot includes a wider collection. 



uncertainties, which in practice means using high 
momentum scales, a penalty of poorer statistics 
must be paid. This effect will become less evident 
if sufficient statistics become available to reduce 
the dominance of the statistical uncertainty, as is 
already the case in some analyses. 
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Figure 15. Comparison of various determinations 
of as with a world average value. 



9. SUMMARY 

In summary, the study of jet distributions at 
HERA, as at the Tevatron, has been a highly 
fruitful area to investigate a number of aspects 
of Quantum Chromo-Dynamics. The theory con- 
tinues to pass all tests, and evaluations of the 
coupling constant as are being performed from 
a number of different perspectives. Accuracy is 
improving with the accumulation of larger data 
samples and further developments in the theoret- 
ical calculations that are available for the experi- 
mentalists to use. However there is an increasing 
need for more NNLO calculations. Nevertheless, 
this is impressive progress and should come to 
a fruition when the full HERA II analyses have 
been completed. 

I should like to thank the organisers of the Ring- 
berg 2008 workshop for all the efforts that they 
put in to ensure an outstandingly enjoyable and 
stimulating meeting. 
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